Joint phenotypes, evolutionary conflict and the fundamental theorem of natural selection

Multiple organisms can sometimes affect a common phenotype. For example, the portion of a leaf eaten by an insect is a joint phenotype of the plant and insect and the amount of food obtained by an offspring can be a joint trait with its mother. Here, I describe the evolution of joint phenotypes in quantitative genetic terms. A joint phenotype for multiple species evolves as the sum of additive genetic variances in each species, weighted by the selection on each species. Selective conflict between the interactants occurs when selection takes opposite signs on the joint phenotype. The mean fitness of a population changes not just through its own genetic variance but also through the genetic variance for its fitness that resides in other species, an update of Fisher\u27s fundamental theorem of natural selection. Some similar results, using inclusive fitness, apply to within-species interactions. The models provide a framework for understanding evolutionary conflicts at all levels

Order-of-magnitude speedup for steady states and traveling waves via Stokes preconditioning in Channelflow and Openpipeflow

Steady states and traveling waves play a fundamental role in understanding hydrodynamic problems. Even when unstable, these states provide the bifurcation-theoretic explanation for the origin of the observed states. In turbulent wall-bounded shear flows, these states have been hypothesized to be saddle points organizing the trajectories within a chaotic attractor. These states must be computed with Newton's method or one of its generalizations, since time-integration cannot converge to unstable equilibria. The bottleneck is the solution of linear systems involving the Jacobian of the Navier-Stokes or Boussinesq equations. Originally such computations were carried out by constructing and directly inverting the Jacobian, but this is unfeasible for the matrices arising from three-dimensional hydrodynamic configurations in large domains. A popular method is to seek states that are invariant under numerical time integration. Surprisingly, equilibria may also be found by seeking flows that are invariant under a single very large Backwards-Euler Forwards-Euler timestep. We show that this method, called Stokes preconditioning, is 10 to 50 times faster at computing steady states in plane Couette flow and traveling waves in pipe flow. Moreover, it can be carried out using Channelflow (by Gibson) and Openpipeflow (by Willis) without any changes to these popular spectral codes. We explain the convergence rate as a function of the integration period and Reynolds number by computing the full spectra of the operators corresponding to the Jacobians of both methods.Comment: in Computational Modelling of Bifurcations and Instabilities in Fluid Dynamics, ed. Alexander Gelfgat (Springer, 2018

Boolean Dynamics with Random Couplings

This paper reviews a class of generic dissipative dynamical systems called N-K models. In these models, the dynamics of N elements, defined as Boolean variables, develop step by step, clocked by a discrete time variable. Each of the N Boolean elements at a given time is given a value which depends upon K elements in the previous time step. We review the work of many authors on the behavior of the models, looking particularly at the structure and lengths of their cycles, the sizes of their basins of attraction, and the flow of information through the systems. In the limit of infinite N, there is a phase transition between a chaotic and an ordered phase, with a critical phase in between. We argue that the behavior of this system depends significantly on the topology of the network connections. If the elements are placed upon a lattice with dimension d, the system shows correlations related to the standard percolation or directed percolation phase transition on such a lattice. On the other hand, a very different behavior is seen in the Kauffman net in which all spins are equally likely to be coupled to a given spin. In this situation, coupling loops are mostly suppressed, and the behavior of the system is much more like that of a mean field theory. We also describe possible applications of the models to, for example, genetic networks, cell differentiation, evolution, democracy in social systems and neural networks.Comment: 69 pages, 16 figures, Submitted to Springer Applied Mathematical Sciences Serie

Identification of hip fracture patients from radiographs using Fourier analysis of the trabecular structure: a cross-sectional study

Peer reviewedPublisher PD

Differential Regulation of the Period Genes in Striatal Regions following Cocaine Exposure

Several studies have suggested that disruptions in circadian rhythms contribute to the pathophysiology of multiple psychiatric diseases, including drug addiction. In fact, a number of the genes involved in the regulation of circadian rhythms are also involved in modulating the reward value for drugs of abuse, like cocaine. Thus, we wanted to determine the effects of chronic cocaine on the expression of several circadian genes in the Nucleus Accumbens (NAc) and Caudate Putamen (CP), regions of the brain known to be involved in the behavioral responses to drugs of abuse. Moreover, we wanted to explore the mechanism by which these genes are regulated following cocaine exposure. Here we find that after repeated cocaine exposure, expression of the Period (Per) genes and Neuronal PAS Domain Protein 2 (Npas2) are elevated, in a somewhat regionally selective fashion. Moreover, NPAS2 (but not CLOCK (Circadian Locomotor Output Cycles Kaput)) protein binding at Per gene promoters was enhanced following cocaine treatment. Mice lacking a functional Npas2 gene failed to exhibit any induction of Per gene expression after cocaine, suggesting that NPAS2 is necessary for this cocaine-induced regulation. Examination of Per gene and Npas2 expression over twenty-four hours identified changes in diurnal rhythmicity of these genes following chronic cocaine, which were regionally specific. Taken together, these studies point to selective disruptions in Per gene rhythmicity in striatial regions following chronic cocaine treatment, which are mediated primarily by NPAS2. © 2013 Falcon et al

Virtual environment navigation with look-around mode to explore new real spaces by people who are blind

Background. This paper examines the ability of people who are blind to construct a mental map and perform orientation tasks in real space by using Nintendo Wii technologies to explore virtual environments. The participant explores new spaces through haptic and auditory feedback triggered by pointing or walking in the virtual environments and later constructs a mental map, which can be used to navigate in real space. Methods. The study included 10 participants who were congenitally or adventitiously blind, divided into experimental and control groups. The research was implemented by using virtual environments exploration and orientation tasks in real spaces, using both qualitative and quantitative methods in its methodology. Results. The results show that the mode of exploration afforded to the experimental group is radically new in orientation and mobility training; as a result 60% of the experimental participants constructed mental maps that were based on map model, compared to only 30% of the control group participants. Conclusion. Using technology that enabled them to explore and to collect spatial information in a way that does not exist in real space influenced the ability of the experimental group to construct a mental map based on the map model

Enrichment analysis of Alu elements with different spatial chromatin proximity in the human genome

Transposable elements (TEs) have no longer been totally considered as “junk DNA” for quite a time since the continual discoveries of their multifunctional roles in eukaryote genomes. As one of the most important and abundant TEs that still active in human genome, Alu, a SINE family, has demonstrated its indispensable regulatory functions at sequence level, but its spatial roles are still unclear. Technologies based on 3C(chromosomeconformation capture) have revealed the mysterious three-dimensional structure of chromatin, and make it possible to study the distal chromatin interaction in the genome. To find the role TE playing in distal regulation in human genome, we compiled the new released Hi-C data, TE annotation, histone marker annotations, and the genome-wide methylation data to operate correlation analysis, and found that the density of Alu elements showed a strong positive correlation with the level of chromatin interactions (hESC: r=0.9, P<2.2×1016; IMR90 fibroblasts: r = 0.94, P < 2.2 × 1016) and also have a significant positive correlation withsomeremote functional DNA elements like enhancers and promoters (Enhancer: hESC: r=0.997, P=2.3×10−4; IMR90: r=0.934, P=2×10−2; Promoter: hESC: r = 0.995, P = 3.8 × 10−4; IMR90: r = 0.996, P = 3.2 × 10−4). Further investigation involving GC content and methylation status showed the GC content of Alu covered sequences shared a similar pattern with that of the overall sequence, suggesting that Alu elements also function as the GC nucleotide and CpG site provider. In all, our results suggest that the Alu elements may act as an alternative parameter to evaluate the Hi-C data, which is confirmed by the correlation analysis of Alu elements and histone markers. Moreover, the GC-rich Alu sequence can bring high GC content and methylation flexibility to the regions with more distal chromatin contact, regulating the transcription of tissue-specific genes

Detecting Linkage between a Trait and a Marker in a Random Mating Population without Pedigree Record

Modern linkage-based approaches employing extended pedigrees are becoming powerful tools for localizing complex quantitative trait loci. For these linkage mapping methods, it is necessary to reconstruct extended pedigrees which include living individuals, using extensive pedigree records. Unfortunately, such records are not always easy to obtain and application of the linkage-based approaches has been restricted. Within a finite population under random mating, latent inbreeding rather than non-random inbreeding by consanguineous marriages is expected to occur and is attributable to coalescence in a finite population. Interestingly, it has been revealed that significant random inbreeding exists even in general human populations. Random inbreeding should be used to detect the hidden coancestry between individuals for a particular chromosomal position and it could also have application in linkage mapping methods. Here we present a novel method, named finite population based linkage mapping (FPL) method, to detect linkage between a quantitative trait and a marker via random inbreeding in a finite population without pedigree records. We show how to estimate coancestry for a chromosomal position between individuals by using multipoint Bayesian estimation. Subsequently, we describe the FPL method for detecting linkage via interval mapping method using a nonparametric test. We show that the FPL method does work via simulated data. For a random sample from a finite population, the FPL method is more powerful than a standard pedigree-based linkage mapping method with using genotypes of all parents of the sample. In addition, the FPL method was demonstrated by actual microsatellite genotype data of 750 Japanese individuals that are unrelated according to pedigree records to map a known Psoriasis susceptible locus. For samples without pedigree records, it was suggested that the FPL method require limited number of individuals, therefore would be better than other methods using thousands of individuals

Design principles for riboswitch function

Scientific and technological advances that enable the tuning of integrated regulatory components to match network and system requirements are critical to reliably control the function of biological systems. RNA provides a promising building block for the construction of tunable regulatory components based on its rich regulatory capacity and our current understanding of the sequence–function relationship. One prominent example of RNA-based regulatory components is riboswitches, genetic elements that mediate ligand control of gene expression through diverse regulatory mechanisms. While characterization of natural and synthetic riboswitches has revealed that riboswitch function can be modulated through sequence alteration, no quantitative frameworks exist to investigate or guide riboswitch tuning. Here, we combined mathematical modeling and experimental approaches to investigate the relationship between riboswitch function and performance. Model results demonstrated that the competition between reversible and irreversible rate constants dictates performance for different regulatory mechanisms. We also found that practical system restrictions, such as an upper limit on ligand concentration, can significantly alter the requirements for riboswitch performance, necessitating alternative tuning strategies. Previous experimental data for natural and synthetic riboswitches as well as experiments conducted in this work support model predictions. From our results, we developed a set of general design principles for synthetic riboswitches. Our results also provide a foundation from which to investigate how natural riboswitches are tuned to meet systems-level regulatory demands